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Abstract. The multiple scattering of photons in a hot, resonant, atomic vapor is investigated and shown to 
exhibit a Levy Flight-like behavior. Monte Carlo simulations give insights into the frequency redistribution 
process that originates the long steps characteristic of this class of random walk phenomena. 

PACS. 02.50.Ey Stochastic processes - 32.50.+d Fluorescence, phosphorescence (including quenching) 



o 

O 

>V 
43 

a* 



> 
o 

cn 
q 

o 
o 



X 



1 Introduction 

Random walk processes of particles, such as the Brown- 
ian motion, are usually described as Gaussian stochastic 
processes, where the distribution P(x) of the size x of the 
random steps has finite first and second moment, so that, 
according to the Central Limit Theorem (in the limit of 
a large number of steps), the total displacement, r, dis- 
tribution converges towards a Gaussian and the process 
is said to be normal or diffusive. Properties such as the 
square root dependence of the total displacement with 
time characterize this so-called normal diffusion (7=1 
in Eq.[Tj): 



< r 2 >= Dt'<. 



(1) 



The 1990's saw the emergence of a new class of sys- 
tems which do not follow the Central Limit Theorem, and 
whose dynamics is dominated by large but rare random 
events. Systems exhibiting long-tailed step distributions 
with diverging second moment give rise to superdiffusion, 
with 7 > 1 in Eq. [T] If this step size distribution asymp- 
totically scales like a power-law P(x) oc x~ a with a < 3, 
the total displacement profile converges towards a Levy 
distribution after a large number of steps [T]. Note that, if 
a < 2, even the first moment of the step size distribution 
diverges, i.e. the mean-free-path between two scattering 
events cannot be defined. 

The number of physical phenomena shown to follow 
Levy instead of Gaussian statistics has been steadily grow- 
ing in the last 20 years |2I3I4| . Recently, Levy processes 
involving light propagation have been demonstrated in an 
optical material engineered to induce photon superdiffu- 
sive behavior [5]. In such a system, long steps take their 
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origin in inhomogeneities of spatial order artificially intro- 
duced in the material. As we shall show, another, much 
simpler, photonic system that exhibits Levy flights, is a 
hot atomic vapor. In that case however, spatial Levy flights 
of resonant light take their origin in spectral inhomo- 
geneities. 

In this paper, we present our study of Levy flights 
of photons in atomic vapors. Contrary to most previous 
work, we measure directly the microscopic ingredient of 
the superdiffusive behavior, i.e. the single step-size distri- 
bution. In the next section, we present the physical origin 
of the Levy nights as well as our experimental methods 
and measurements. Then, in section [3j we detail our nu- 
merical simulations of such a system. 



2 Photon diffusion in atomic vapors 
2.1 Origin of long steps 

Radiation trapping in an atomic vapor is a well known 
phenomenon [6 7[Sj, whereby atoms in a vapor absorb res- 
onant radiation from an incident source and re-emit pho- 
tons that can be re-absorbed in the vapor, and so on. 
This photon multiple scattering happens in atomic media 
ranging from laboratory atomic vapors in optical cells to 
interstellar space. Under simple hypotheses, one can show 
that the distribution for the length covered by a photon 
between an emission and an absorption (jump size dis- 
tribution) asymptotically decays like a power law of di- 
verging second moment, characterizing resonant photons 
trajectories in atomic vapors as Levy flights j9]. The diffu- 
sion of photons exhibits an anomalous behavior through 
the inhomogeneity introduced in the scattering process by 
the frequency dependence of the absorption coefficient and 
the frequency redistribution due to the atoms velocity. 
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In a dilute atomic vapor of two level atoms where the 
Doppler broadening of the atomic lineshape is negligible in 
comparison with the natural width, as with laser-cooled 
atomic vapors [10111112113] . the scattering is quasi elas- 
tic, i.e. the frequency of the re-emitted photon is almost 
the same as the frequency of the absorbed one, so that 
frequency redistribution can be neglected. For light at a 
frequency the jump size distribution then writes 



dT(x, u) 
dx 



k a {v)e 



-k a {v)x 



(2) 



where k a {v) is the absorption profile and T{x, v) the prob- 
ability of a photon of frequency v going through a distance 
x without being absorbed, 



T(x,u) = e-M^ 



(3) 



In the general case where inelastic scattering cannot be 
discarded, as in normal laboratory temperature condition, 
however, the jump size distribution function has to be 
computed by averaging over the spectrum of emitted light 
k e : 

^ k e {v)k a {v)e- k ^ x dv . (4) 



P(x) = [ 
Jo 



In presence of Doppler broadening alone, the absorp- 
tion profile is gaussian 



k D (v) = Na 



r 



-xH 2 



while it is lorentzian for natural broadening 



1 



k L (v)=Ncr 

i + IM" - w)ir\ 

and Voigt if the two mechanisms are mixed, 



(5) 



(6) 



r a 

4 7T 



exp(-r 



ttvqVo J.oo a 2 + (x(f) - y) 



2 dy 



Here, N is the atomic density, <r$ the on-resonance atomic 
cross section for atoms at rest, vq the resonance frequency, 
and r the natural width of the atomic transition. The 



2.2 Experimental methods and measurements 

While the measurement of the single step size distribution 
of a photon in large, hot atomic vapors such as in stars 
is not feasible, a specific multicell arrangement allowed 
for such an experimental observation [14] . The experimen- 
tal arrangement devised to observe individual scattering 
events in a hot vapor consists of an optical cell filled with 
a vapor of Rb atoms at a controlled temperature/density, 
from now on called the observation cell (part (a) of Figure 
[1} . These atoms are excited by a collimated beam of pho- 
tons, so as to define a direction of observation (the ballis- 
tic one). The atoms scatter the incident photons in all the 
directions. Photons detected by a CCD device show the 
position in the cell where they have been last scattered. A 
narrow strip (2 mm) parallel to the ballistic direction is se- 
lected for measurements. Assuming that photons scattered 
from this strip to outside the cell have undergone a sin- 
gle scattering event in the observation cell, the measured 
fluorescence intensity is directly proportional to the step- 
size distribution P{x). In fact, insuring a single scattering 



most probable velocity at temperature T is vq = ^j2ksT jm 
with ks the Boltzman constant, c the speed of light and 
m the atomic mass. x{ u ) = ^ u ~^° is a reduced frequency. 
Finally, to define the Voigt profile, we have introduced the 
Voigt parameter a = / n c/(47ruoi'o). The determination of 
the emission profiles, however, is not straightforward, as 
they depend on the light initial spectra, but Monte-Carlo 
simulations show that, in the multiple scattering regime, 
they are very close to the absorption lines (see Section [3]). 

For Gaussian emission and absorption spectra k a {v) = 
k e {y) = fcrj(^), the expected power of the step length dis- 
tribution asymptotic decay is a = 2, while for Lorentzian 
(Cauchy) spectra k a {u) = k e {v) = fci(^), a = 1.5. As 
Voigt spectra ky{v) show Lorentzian-like frequency wings, 
they also lead to a decay with a = 1.5 [5]. 
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Fig. 1. Experimental setup for observing the spatial profile 
of fluorescence in a hot atomic vapor, i.e. the step-size dis- 
tribution function, (a) Observation cell and imaging system. 
This part is common to the three performed measurements, 
(b) The incident photons are directly provided by a monochro- 
matic laser, for calibration, (c) The incident photons have been 
scattered once, at 90°, in the low-optical-thickness source-cell, 
(d) The incident photons have been first scattered several times 
(~ 4) in a first, optically-thick, cell before being scattered once 
more in the second cell. 
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regime in the observation cell is quite incompatible with 
having a sufficient dynamics of P(x) so that higher order 
scattering inevitably happens. As multiple scattering cor- 
rections might also mimic a power law intensity along Ox, 
it is crucial to correct for this effect. This can be achieved 
by subtracting a narrow strip measured slightly off-axis, 
as those photons are due to multiple scattering only (14j . 

A reference measurement is first carried out with a 
monochromatic laser beam (frequency i>l) as the incident 
photon source (See FigureQJb)). The step-size distribution 
P(x) for this specific configuration is given by Eq. [5] with 
v = Vl , 

P{x) oc e - fc (^)*. (8) 

The measurement of this exponential decay enables the 
calibration of the mean-free-path £(vl) = l/fc(^) at res- 
onance, and consequently the atomic density n, which 
can be varied through the temperature of the Rb reser- 
voir in the cell. The temperature is varied between 20 
and ~ 47° C, corresponding to an atomic density between 
9 x 10 9 and 2 x 10 11 atoms/cm 3 and a subsequent resonant 
mean free path between 50 and 5 mm. 

A second experimental measurement of the jump size 
distribution P(x) is performed by using a two cell config- 
uration (Fig. (He)). Laser light is sent on a first rubidium 
cell of very low atomic density, so that photons are scat- 
tered at most once in this source cell. Light with a spectra 
broadened by Doppler effect is emitted in all directions in 
space, but a collimated beam orthogonal to the initial laser 
propagation axis is selected out of it by two diaphragms 
and sent towards the observation cell. The knowledge of 
the position of the first scattering event, in the source cell, 
and the second one, in the observation cell, gives us access 
to the jump size distribution function P(x) (Fig. We 
find that, for x large enough compared to the mean- free- 
path of resonant light, P{x) scales like a power law l/x a , 
with a = 2.4f ±0.12. 

Although the obtained distribution has clearly a di- 
verging second moment, this measurement alone does not 
prove that light transport in the multiple scattering regime 
is superdiffusive, because the spectrum of light scattered 
in the first cell, and hence the jump size distribution for 
the first step, depends on the laser frequency. The fre- 
quency redistribution is only partial [15] . To characterize 
the multiple scattering regime, we need to measure the 
jump size distribution function after several steps, once 
photons have lost memory of the initial laser frequency, so 
that the emission spectra has converged towards a limit 
spectra (see section [3] for a numerical demonstration of 
this convergence). To achieve this, we use a three cell con- 
figuration (Fig. HJd)). Laser photons are sent in a first 
rubidium cell of high atomic density, where they are scat- 
tered several times (~ 4), before reaching a second source 
cell, where they undergo one more scattering event with a 
well-known position before being sent to the observation 
cell. We thus measure the jump size distribution func- 
tion in the multiple scattering regime. Once again, we 
find a power law asymptotic behavior P(x) ~ l/s", with 
a = 2. 09 ±0.15, which has an infinite variance and is there- 
fore characteristic of a superdiffusive regime. One should 
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Fig. 2. Experimental measurements of the step size distribu- 
tion function (log-log scale) in the configuration of Figure [Ha) 
and (b). With an incident laser beam (green crosses), the ob- 
served decay is well fitted by an exponential (red dotted line) 
as expected from Beer-Lambert law. In the two-cell configura- 
tion, the data (continuous blue line) are well fitted by a power 
law (black dotted line) P(x) = l/x a , with a = 2.41 ± 0.12. 



notice, however, that as a > 2 the mean free path remains 
finite. 



3 Monte Carlo simulations 

A Monte Carlo (MC) code simulates the spectral and spa- 
tial properties of individual photons as they propagate in 
an atomic vapor at temperature T, in an optical cell hav- 
ing the dimensions of the experimental observation cell. 
Monte Carlo simulations may in general help in the sig- 
nal deconvolution or in choosing suitable experimental pa- 
rameters. In this work, it proved indeed useful to validate 
the data analysis and for a deeper understanding of the 
third experimental configuration, in the multiple scatter- 
ing regime (Fig. [Ijd)). 



3.1 Description 

The code consists essentially in an elementary procedure, 
corresponding to a scattering order, which is repeated as 
long as the simulated photon is inside the optical cell. 
When it gets out, it is either discarded or detected, de- 
pending on the direction it takes. If detected, its last scat- 
tering position in the cell is recorded in a file, together with 
the order of the last scattering process. The elementary 
procedure consists in taking one photon, of frequency v\ 
and wavenumber ki , and to determine the position where 
it will undergo a scattering process. The distance from 
this position to the initial one is frequency-dependent. The 
probability T{x, v) that a photon of frequency v travel a 
distance x without being absorbed in the vapor is ran- 
domly attributed a value between and 1 and the step 
length x is then determined as 



InT 



(9) 
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with k(v) the absorption spectrum. A new wavenumber 
k.2 (direction) is then drawn for the scattered photon and 
its frequency v 2 is calculated from 



The simulated first-step distribution is shown in Figure 
[3] to reproduce the observed one. 



V-2 



v\ + v • (k 2 - ki), 



(10) 3.2 Multiple scattering regime 



where v is the velocity of the scattering atom. The fre- 
quency will then determine the step length to the next 
scattering process. 

The photons incident on the optical cell are initially 
homogeneously distributed over a disk whose diameter is 
the one of the second diaphragm of the setup shown in 
Figure l(c,d). Their direction is parallel to the longitudi- 
nal axis of the observation cell and their frequency can be 
drawn either from a laser's lorentzian lineshape (calibra- 
tion experiment, see text) or from a Voigt distribution at 
the temperature of the source cell. 

From the equation (fTDf . we can see that the relevant 
components of both the wavenumbers and of the atom 
velocity are the one parallel and the one in the plane per- 
pendicular to the incident photon's direction. Regarding 
the direction of the scattered photon, an isotropic angular 
distribution follows from isoprobabilistically drawing the 
value of the polar angle's cosine between —1 and +1 and 
the value of the azimuthal angle between and 2ir. 

Finally, the velocity component v p of the scattering 
atom in the plane perpendicular to the direction of the 
incident photon can take any value, with a probability law 
given by the Gaussian Maxwell-Boltzmann distribution of 
standard deviation <x„ = ^Jk B T/m: G(v p ) cx e~ < - Vp /^ 2 ' T ^ . 
The velocity component Vk of the scattering atom in the 
direction of the incident photon is drawn from a prob- 
ability law p(ffc) = L(v,Vk) x G(wfc), convolution of the 
atomic Lorentzian lineshape L(y, Vk) cx 1/(1 +4tt(is — vq — 
Vk / X) 2 / r 2 ) and of the Gaussian Maxwell-Boltzmann ve- 
locity distribution G(vk), where vq, A and r are respec- 
tively the frequency, the wavelength and the linewidth of 
the atomic transition considered. 



1000 




x (cm) 

Fig. 3. First-step size distribution function (log-log scale) in 
the configuration of Figure [T] Black solid line: Experimental 
measurements, as in Figure [21 Grey line: Monte Carlo simula- 
tion. 



To characterize the multiple scattering regime from the 
measurement of a single step-size distribution in the mid- 
dle of the scattering sequence, it is necessary that this 
distribution have reached a stationnary state and do not 
depend any more of the photon frequency at the previous 
step. To know how many scattering events are required 
to fulfill this condition, we have used the MC simulations 
to follow the evolution of the light spectral distribution, 
from the monochromatic laser beam (lorentzian shape of 
width sa 1 MHz) to the asymptotical Voigt distribution 
(Doppler width for the preparation cell temperature) . The 
half width at half maximum (HWHM) of the spectra is 
shown as a function of the diffusion order in Figure Ufa). 
We can see that in the conditions of the experiment (300 
K) at least 5 diffusions are necessary for assuring a Voigt- 



(a) 




4 6 
Scattering events 



(b) 




Scattering events 

Fig. 4. Evolution of (a) the HWHM of the photons spectral 
distribution and (b) the power-law exponent a, with the num- 
ber of scattering events in a Rb vapor at 300K. The original 
spectral distribution (photons absorbed in the zeroth order dif- 
fusion) is the one of a resonant, monochromatic laser. 
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broadened distribution of the photons incident on the ob- 
servation cell and fulfill the Complete Frequency Redis- 
tribution condition. The second setup used in the exper- 
iments (see Figure Hid)) takes this requirement into ac- 
count. The exponent of the power-law of the correspond- 
ing step-distribution is shown in Figure Q|b). It decreases 
from a ss 2.4 to ss 2.05 as the diffusion order increases 
from 1 to ~ 5. 

In our experimental conditions, the Doppler width is 
typically two orders of magnitude larger than the natural 
linewidth (Voigt parameter a « 10 -2 ). In such a regime, 
the lorcntzian nature of the natural broadening is expected 
to dominate the spectral dynamics (asymptotical Lorentz 
power-law exponent a = 1.5) for optical densities na x 
larger than 10 4 [S], i.e. for jump sizes larger than 1 m in 
our experimental conditions (cto ~ 1.25 x 10 -13 m 2 , n ~ 
5 x 10 16 m~ 3 ). Below 10 cm, the spatial range observed 
in our experiment, the spectral behavior (assuming CFR) 
is purely Doppler-driven, with an exponent a = 2 for the 
step distribution power law. 

4 Conclusion 

We have evidenced a Levy Flights behavior of photons 
propagation in a hot resonant atomic vapor by measur- 
ing the single step distribution in this multiple scattering 
process. Monte Carlo simulations help deconvoluating the 
effects of multiple scattering and analyzing the frequency 
redistribution process in a Doppler-broadened vapor ex- 
cited by a monochromatic light source. 
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